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ABSTRACT 

The disk galaxy UGC8802 has high neutral gas content and a flat profile of star formation 
rate compared to other disk galaxies with similar stellar mass. It also shows a steep metallicity 
gradient. We construct a chemical evolution model to explore its growth history by assuming 
its disk grows gradually from continuous gas infall, which is shaped by a free parameter - the 
infall-peak time. By adopting the recently observed molecular surface density related star for- 
mation law, we show that a late infall-peak time can naturally explain the observed high neutral 
gas content, while an inside-out disk formation scenario can fairly reproduce the steep oxygen 
abundance gradient. Our results show that most of the observed features of UGC8802 can be 
well reproduced by simply 'turning the knob' on gas inflow with one single parameter, which 
implies that the observed properties of gas-rich galaxies could also be modelled in a similar way. 

Subject headings: Galaxies:evolution - Galaxies:photometry - Galaxies: stellar content 



1. Introduction 

Understanding the content and distribution of 
cold gas in galaxies is an important step to un- 
derstand the formation and evolution of galax- 
ies. UGC8802 is an interesting target selected 
from the GALEX Arecibo SPS S Survey (GASS) 



(iGarcia-Appadoo et al. 20091 : Catinella et al 



20101 ). One interesting aspect is its cold gas con- 
tent. UGC8802 is a disk galaxy with stellar mass 
M* « 2 x 10 10 M Q , while its HI mass is estimated 
to be as high as 2.1 x 10 10 M© (jSpringob et al. 



20051 ) and its molecular ga s mass is only abou t 



one-tenth of its HI mass (jMoran et al. I l2010h 



Compared to other disk galaxies in this stel- 
lar mass range, the high neutral gas content of 



UGC8802 is uncommo n dGiovanelli et al.l 120071 



Garcia- Appadoo et al. I l2009t iCatinella et al 



20101 ). Why UGC8802 has such high neutral gas 
fraction and how UGC 8802 acquires its col d gas 



are still open questions ( Moran et al. 2010l) 



Another interesting aspect of UGC8802 is 
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i ts flat profile of th e star formation rate (SFR) 



(jMoran et al. Il2010l ). The SFR surface densities 



of 'normal' disk galaxies typically decrease (some- 
times exponentially) from the center to the outer 
regi ons of the disk, such as the Milky Way, M31, 
etc (iBieiel et al. ll2008tlFu et al. lliooalYin et al. 
l2009h . Besides the fiat SFR profile. iMoran et al. 
<|2010h found that UGC8802 shows a steeper ra- 
dial oxygen gradient than the average gradient of 
the GASS sample. 

Therefore, it is interesting to explain these ob- 
served results under the frame of the current mod- 
els of galaxy formation and evolution. For individ- 
ual disk galaxy, parameterized modeling is proven 

thei r galactic formation 



to be fruitful in exploring thei 
and evolution ( Tinslcv Ill980l 



Chang et al 



Boissier fc Prantzos ll2000tlC hiappin i et al 



1999: 



2001 



Colavitti et al. 1 120091: lYin et aL Fl2009h . In this 
letter, we construct such a parameterized model 
to investigate the radial-dependent star formation 
history (SFH) of UGC8802 and so as to explore 
its growth history. 

2. The model 



Our model is based on that of iChang et al 



(<2010h . Here we emphasize main ingredients of the 



1 



model, especially some revisions. We adopt the 
standard cold dark matter (CDM) cosmology with 
n M = 0.3, n A = 0.7 and H = 70kms- 1 Mpc- 1 . 

We assume that the galactic disk is sheet-like 
and composed by a set of independent rings, each 
500pc wide. The disk originates through continu- 
ous gas-infall from the dark halo. We assume that 
the disk begins to form at 12.5 Gyr ago, which 
roughly corresponds to z ~ 6 under the standard 
cosmology. At given radius r, the gas infall rate 
/i n (r, t) (in units of M f ^pc~ 2 Gyr~ 1 ) is a s sume d to 
be Gaussian in time (jChang et al.lll999ll2010h : 



/inM) 



A(r) 



-[t-t p (r)] 2 /2a 2 



2na 



(1) 



where t p is the infall-peak time and a is the scat- 
ter. We adopt t p to be free and set a = 3Gyr 
since the model results are insensitive to v arying 
cr in the 2 ~ 4Gyr range (|Chang et al.ll2010h . The 
A(r) are a set of separate quantities normalized 
by the stellar mass surface density of present-day 
E*(r, t g ), where t g is the cosmic age and set to be 
t g = 13. 5 Gyr according to the adopted cosmology. 
The stellar mass of UGC8802 and its disk scale- 
length at the present time are adopted to be M* 



2 x 1O 1O M and r d = 5.8kpc (jMoran et al. Il2010l) 



We assume that the galaxy UGC8802 at present is 
a pure-disk system and t g ) follows an expo- 
nential profile, and so that the central stellar mass 
surface density is given as E*(0,i g ) = M*/{2-Kr\). 

The star formation (SF) law is one of the key 
ingredients of ou r disk galaxy formation model. 
Kennicutt I (|l998h found a power-law relationship 
between the galaxy-averaged SFR surface density 
and the galaxy-averaged total gas surface density, 
which is often called as the classical Kennicutt- 
Schmidt SF law and has been widely used in the 
chemical evolution models of the disk galaxy for- 
mation. Later, extended analyses to kpc-scale 
regions within several Local Group galaxies show 
that the SFR surface density correlates stronger 
with the surface density of mole cular hydrogen 



than that of atomic hydrogen (I Wong k, Blitz 



2002t iBigiel et al. II2008I: iLerov et al. 1120081 ) 



In this study, we adopt a recently observed 
molecular surface density related SF law (here- 
after, it is termed as the Sn 2 -based SF law), that 
is, the SFR surface density *(r, t) is linearly pro- 
portional to the molecular hydrogen surface den- 



sity £ H2 (r,i): 



*(r,t) = S H2 (r,t)/t de p, 



(2) 



where tc\ PV is the m olec ular gas depl e tion t ime 



Bigiel et all (|2008h and ILerov et all (|2008l ) de 



rived a constant molecular gas d epletion time as 
tdm = 2Gyr in their sample, while lSaintonge et al. 
(|201ll ) found the non-universality of idep in galax- 



ies in the CO Legacy Database for the GASS 
(COLD GASS) sample. ISaintonge et~ai~l (l201ll) 



is 



found that the strongest dependence of idep 
on the stellar mass and the mean idep is pa- 
rameterized as log(£dep/yr) = 0.36(logMH,/M Q — 
10.7) + 9.3. According to this relation, we adopt 
t dcp = 0.77Gyr for UGC8802. 

Regarding th e molecular-to-atomic mass frac- 
tion R mol (r,t), iBlitz fc Rosolowskv I (|2006l ) ob- 
tained an empirical fit relation: 

R mol (r,t) = EH a (r,f)/Em(r,t) = {P h /Ph,o)\ 

(3) 

where is the mid-plane pressure, Ph.o an d 7 
are fitting constants and we adopt Ph.o/k = 1.7 x 
10 4 cur 3 K and 7 = 0.8 (|Lerov et al. II2008I) . The 
pressure is estimated to include gas self-gravity 
(|Elmegreenlll989HLerov et al. 1120081) : 



Pi, 



(4) 



where £„ 



gas is the total gas surface density, a g and 
<j* iZ are the vertical velocity dispe rsion of gas and 
stars . We adopt a g = llkms -1 (jOstriker et al. 
2010h. while r 



is es timated as [see appendix B 
of ILerov et al. I (|2008h for details]: 



(5) 



A similar EH 2 -based SF law has been suc- 
cessfully adopted i n the semi-analyt i cal model of 
galaxy formation by lFu et al~ ( 201d( ). Kang et al. 
(|2012l ) also use the Eh 2 -based SF law to explore 
the chemical evolution and SFH of M33 and find 
that, comparing to the model adopting the classi- 
cal Kennicutt-Schmidt SF law, the model adopt- 
ing the Eh 2 -based SF law predicts steeper color 
and matellicity gradients, which are more consis- 
tent with the observations. 

Other components of the model, such as the ini- 
tial mass faction (IMF), chemical evolution, etc., 
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are the same as that of Chang et al. ( 2010h . We 
also consider the contribution of gas outflow pro- 
cess by assuming that the gas outflow rate is pro- 
portional to the SFR and the coefficient is set to 
be 6 ut = 0.007 according to the mass-dependent 
model of Chang et al. ( 20f0t ). We point out that 
our final results are insensitive to the variation of 
6 out since UGC8802 is a massive galaxy and the 
gas outflow process is not a dominant process of 
its chemical evolution. 

We emphasize that our model onl y has one free 
param eter, the infall-peak time t p 



(|20f0h . t p is assumed to be a function of galac- 
tic stellar mass and do not vary with radius. 
In the case of UGC8802, we assume that t p is 
a function of radius since we will further con- 
sider the radial profiles of the UGC8802 rather 
than the global pr operties of galaxies as that in 
Chang et al. ( 2010l ). Since the formalizing of t p (r) 
regulates the shape of SFHs along the disk, we aim 
to get the constraints on t p (r) from the observed 
properties of UGC8802 and then get insights on 
its growth history. 

3. Observations versus model predictions 

We present our results step by step. Firstly, 
we summarize main observational properties of 
UGC8802, and then explore the influence of the 
free parameter on our results. Finally, we present 
a viable model to discuss the radial-dependent 
SFH of UGC8802. 

The observed data of the rad ial profiles of 



UGC8 802 are adopted mainly from lMoran et al. 
( 2010h and plotted in Fig. 1. Since our model 



only has the galactic-center distance and can- 
not distinguish the negative position from the 
positive one, we plot the observed data of the 
negative and positive positions as open triangles 
and filled squares in Fig. 1, respectively. The 
NUV-r color profile of UGC8802 is measured us 



ing S PSS and GALEX photometry (jMoran et al 
2010), but the data has not been corrected for 
dust-attenuation. The observed data of oxygen 
abundance log(0/H) + 12, the SFR surface den- 
sity and D„(4000) are obtai ned from the long-sli t 



spectroscopy of UGC8802 (jMoran et al. I l2010h . 
The data SFR/ < SFR > are estimated from 
the observed SFRs and the stellar surface density 
£»(r, t g ) by assuming (1 — R) fraction of the total 



formed stellar mass are locked in the stellar mass 
in the present day, where R is the return fraction 
and we set R = 0.3 according to the adopted IMF. 

The comparison between model predictions of 
the radial profiles and the observations are pre- 
sented in Fig.l. The dash and dash-dot lines 
are model results of two limiting cases of t p = 
O.lGyr and t p = 15Gyr, respectively The case 
tp = O.lGyr (dash lines) corresponds to a time- 
decreasing gas-infall process that most of the gas 
has been accreted in the early stage of its history, 
while that of t p = 15Gyr (dash-dot lines) repre- 
sents a time-increasing gas accretion that there is 
still a large fraction of cold gas-infall at the present 
time. Fig. 1 shows that the model predictions 
are very sensitive to the adopted t p . The model 
adopting the later infall-peak time (dash-dot lines) 
predicts the higher surface density of H2 and HI, 
the higher SFR and specific SFR, the lower gas- 
phase metallicity, the older mean stellar age and 
the bluer color. These results are straightforward, 
the later t p corresponds to a more recent gas infall 
and then results in a younger stellar population 
and a higher fraction of cold gas in present. 

Fig. 1 also shows that the area between the 
dash and dash-dot lines almost covers the loca- 
tions of the observed data. There is a trend 
that the observed data at the inner region ap- 
proaches the dash lines (the early t p ) and the 
data at the outer regions approaches the dash- 
dot lines (the late t p ). After a set of calcula- 
tions and comparisons, we select a viable model 
by adopting an inside-out formation scenario 
t p (r)/Gyr = l.hr/r^ + 5.0 and plot its model re- 
sults as solid lines in Fig. 1. This model predicts 
Af H2 = 3.7 x 1O 9 M and M m = 1.5 x 10 10 M Q , 
which in good agreement with the observed val- 
ues [M HT = 2.1 x 10 10 M m , M H „ = 1.45 x 1Q 9 M 



(Mo ran et al. 



20101) ]. considering the observed 
uncertainties. Our model-predicted NUV — r col- 
ors are systematically bluer than the observed 
ones, which is also reasonable since the observed 
NUV — r color has not been corrected for dust- 
attenuation and our model predictions are dust 
free. 

The good consistence between our model pre- 
dictions and the observations implies that the con- 
tinue gas-infall model is also viable for the 'pe- 
culiar' UGC8802, which has steeper metallicity 
gradient, flatter SFR profile and extraordinary 
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Fig. 1. — Comparisons of the radial profiles be- 
tween the observations and the model predictions. 
The dash and dash-dot lines are model results 
adopting t p = O.fGyr and t p = 15Gyr, respec- 
tively. The solid lines plot the predictions of the 
viable model, which adopts i p (r)/Gyr = 1.5r/ra + 
5. The observed data are described in details in 
the text. 

higher neutral gas fraction than normal disk galax- 
ies. The key ingredient of this successful model 
is that the infall-pcak time t p (r) increases from 
5Gyr in the innermost disk to 12.5Gyr at the out 
region around ~ 5rd- That means the formation 
of UGC8802 is quite late and the disk growing is 
still ongoing. In fact, this scenario is the well- 
known idea that the disk forms insidc-out and has 



and evolution of disk 


galaxies (Chane et al. 


1999; 


Boissier & Prantzos 


2000; ChiaDDini et al. 


2001; 


Colavitti et al. 2009 


; Fu et al. 1 2009; Yin et al. 


200alWanff et al. 1 201lh. 



The above results have shown that the inside- 
out growth scenario can nicely reproduce the main 
observed properties of UGC8802. In Fig. 2, we 
show further predictions of the growth history of 
the UGC8802 of this established model. The up- 
per panel of Fig. 2 shows the model predicted time 
evolution of the half-size stellar mass R50,* (solid 
line), which is defined as the radius where con- 
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Fig. 2.— The growth history of UGC8802 pre- 
dicted by the viable model. In the upper panel, 
the solid line shows the model predicted time evo- 
lution of the half mass size R50,*, while the dash 
line shows the observed trend of size evolution of 
massive galaxies. The bottom panel shows the 
normalized evolution history of different baryonic 
components, including the stellar mass (the solid 
line), the atomic hydrogen mass A/hi (the dash 
line) and the molecular hydrogen mass Mh 2 (the 
dot line). 

tains half of the total stellar mass. As we can see, 
our model predicts that the half-size of UGC8802 
increases almost linearly with time, which is also 
globally consistent with the observations of the av- 
erage size evolution trend of t he massive galaxies 
\rjz) /rjz = 0) = (1+z)- 11 , (|van der Wei et al. I 



20081 ). dash line]. This size growth is actually an- 
other angle of view of the disk inside-out formation 
scenario. 

In the bottom panel, we plot the model predic- 
tions of the evolution of the different mass com- 
ponents, which includes the stellar mass M* (solid 
line), the atomic hydrogen mass Mhi (dash line) 
and the molecular hydrogen mass Mh 2 (dot line) . 
The quantities are normalized by its value at the 
present-day. The bottom panel of Fig. 2 shows 
that, in the early epoch, with the continue in- 
fall of cold gas (mainly in the form of the atomic 
gas), the molecular gas accumulates faster, while 
the accumulation of stellar mass is lagged behind. 
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Later, the mass of the molecular gas becomes to be 
asymptotic, while the growth rate of stellar mass 
is speeded. Indeed, half of the stellar mass has 
been accumulated during the last 4Gyr. To sum 
up, our results suggest that UGC8802 is a young 
galaxy and still active in star formation processes. 

4. Summary 

In this work, we build a bridge for disk galaxy 
UGC8802 between its observed properties and its 
growth history by constructing a simple chemical 
evolution model. We find that a late infall-peak 
time t p , which corresponds to a recent gas-infall 
process, is necessary to explain the observed high 
atomic and molecular gas content of UGC8802. 
The model adopting an inside-out formation sce- 
nario £ p (r)/Gyr = 1.5r/r<j + 5.0 and the £hv 
based SF law can reproduce the observed radial 
profiles of UGC8802 very well. This consistency 
does encourage us that our simple model is viable 
for such a 'peculiar' galaxy. This result implies 
that there should be no violent star burst process 
(e.g. the major merging) during the growth his- 
tory of UGC8802. According to our model, half of 
the stellar mass of UGC8802 may have been accu- 
mulated during the last 4Gyr and the SF process 
of the outer disk is still active. The physics behind 
the very young history of UGC8802 may related 
to its environment, however, which is beyond the 
scope of this letter. 

The success of our simple model on the case 
of UGC8802 implies that the observed properties 
of the gas-rich galaxies might also be well mod- 
eled by 'tuning the knob ' on gas inflow with a sin- 



gle parameter Ind eed, IWang et al. I ([201 II ) and 



M oran et al. I ([20121 ) find that the color and mctal- 
licty of the outer region of the galaxies strongly 
correlate with its total atomic mass Mm, i.e., the 
higher the Mm, the bluer the color, the lower the 
metallicity. In our story, the high atomic gas frac- 
tion of gas-rich galaxy may be the result of a late 
and smooth gas accretion history, and the bluer 
color and lower metallicity of its outer region are 
embedded in the disk inside-out formation sce- 
nario. We will further explore these questions in 
a forthcoming paper. 
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